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Government.  Neither  the  United  States  Government  nor  any  agency  thereof,  nor  any  of  their 
employees,  makes  any  warranty,  express  or  implied,  or  assumes  any  legal  liability  or 
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recommendation,  or  favoring  by  the  United  States  Government  or  any  agency  thereof.  The  views 
and  opinions  of  authors  expressed  herein  do  not  necessarily  state  or  reflect  those  of  the  United 
States  Government  or  any  agency  thereof. 
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EXECUTIVE  SUMMARY 


This  research  focused  on  obtaining  an  understanding  of  the  nanoscale  physical  mechanisms  that 
govern  interfacial  interactions  and  heat  transfer  in  materials  at  elevated  temperatures.  The  underlying 
strategy  was  to  understand  how  thermal  conduction  at  nanoscale  gets  affected  by  factors  such  as 
the  effect  of  quantum  to  molecular  scale  temperature  related  phase  change,  effect  of  nanoscale 
interfacial  reconstruction,  and  effect  of  thermal  stress  affected  material  transformation.  Issue  is  to 
enhance  thermal  management  characteristics  of  materials  exposed  to  high  heating  rates  and  high 
thermal  gradients  based  on  understanding  gained  in  this  research  as  well  as  to  leverage  this 
understanding  to  resolve  thermal  shock  and  thermal  cycling  affected  failures  in  materials 
exposed  to  temperature  extremes  by  introducing  control  factors  related  to  nanoscale  thermal 
conduction.  Material  systems  of  focus  are  Si  based  ceramics,  particularly  ZrB2-SiC 
nanocomposite  interfaces.  For  verification  purposes  analyses  also  focus  on  Si-Ge  interfaces  and 
nanocomposites,  as  a  lot  of  simulation  data  is  available  on  this  system. 

In  terms  of  modeling  the  required  phenomenon,  classical  analytical  approaches  could  not  be 
used  because  of  small  length  scales  at  material  interfaces.  Classical  molecular  dynamics  or 
Monte  Carlo  methods  based  approach  could  not  be  used  because  classical  interatomic  potential 
used  in  such  approaches  are  not  capable  of  correctly  describing  phase  transformation,  high 
temperature  property  changes,  and  electronic  properties  which  must  be  correctly  described  in 
order  to  properly  understand  such  phenomenon.  There  have  been  significant  advancements  in 
experimental  measurements  of  thermal  conduction.  However,  such  advancements  have  not 
addressed  in-situ  high  temperature  thermal  conduction  as  well  as  thermal  conduction  under 
applied  stress  and/or  strain.  Therefore,  developments  were  made  both  in  terms  of  modeling  and 
experiments  in  this  research  to  model  the  required  material  thermal  behavior.  A  new 
experimental  approach  was  developed  to  perform  in-situ  measurements  of  temperature 
dependent  thermal  conduction  behavior  as  a  function  of  applied  stress/strain.  The  experimental 
developments  were  coupled  with  new  developments  in  non-equilibrium  Green’s  function  based 
ab-initio  calculations  of  nanostructure  thermal  conductivity. 

Using  ab-initio  modeling  it  was  revealed  that  interatomic  separation  in  semiconducting  as 
well  as  high  temperature  materials  in  single  crystalline  state  as  well  as  in  interfacial 
configurations  is  intricately  linked  to  change  in  both  electronic  and  phononic  (lattice)  density  of 
states,  especially  at  high  temperatures.  Such  relation  occurs  to  different  extents  in  metal, 
semiconductor,  and  ceramic  single  crystals.  Ab-initio  simulations  were  used  to  highlight  important 
role  played  by  electronic  thermal  conductivity  in  overall  thermal  conduction  across  interfaces  at  high 
temperatures.  This  is  the  first  time,  quantum  simulations  of  electronic  and  phononic  thermal 
conductivity  across  interfaces  in  a  realistic  material  interface  were  reported.  Ab-initio  analyses  have 
also  revealed  new  insights  regarding  role  of  temperature  and  stress  in  phonon,  electron,  and 
coupled  phonon-electron  affected  heat  conduction  in  model  metal  (Au),  semiconductor  (Si),  and 
ceramic  (SiC)  materials.  It  was  shown  that  electronic  thermal  conductivity  becomes  a  significant 
part  of  a  material’s  thermal  conductivity  at  elevated  temperatures.  Calculations  were  used  to 
understand  unresolved  experimental  findings  regarding  electronic  thermal  conductivity  issues  and 
related  failure  issues  in  high  temperature  materials  on  experiments  performed  by  Fahrenholtz 
group  and  Dave  Marshall’s  group.  Discussions  are  underway  on  how  such  analyses  be  used  to 
aid  directed  processing  of  high  temperature  thermal  shock  resistant  ZrB2  based  composites.  The 
analyses  are  also  underway  to  predict  quantum  mechanistic  ways  to  reduce  thermal  shock  and 
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thermal  cycling  related  failure  in  such  materials.  We  have  also  performed  classical  molecular 
simulations  of  thermal  conduction  in  nanoscale  microstructures  in  order  to  examine  scale  up  of 
the  findings  of  classical  ab-initio  simulations  to  realistic  system  sizes.  We  have  proven  for  the 
first  time  that  materials  with  biomimetic  phase  morphology  have  thermal  conductivity  values 
independent  of  strain.  This  finding  has  strong  implication  for  developing  materials  with  thermal 
properties  independent  of  applied  stress.  We  have  proven  for  the  first  time  that  that  tensile 
straining  and  heat  flow  direction  can  be  used  to  develop  a  thermal  diode  material  from 
superlattice  construction.  In  addition,  our  group  was  the  first  to  show  that  nanostructures  with 
tunable  thermal  diode  properties  could  be  developed  based  on  strain  engineering.  For  this  reason, 
the  relevant  publication  was  featured  in  the  Virtual  J.  of  Nanoscale  Science  and  Technology  (A 
collection  of  significant  advances  in  nanotechnology). 

Experimental  analyses  have  revealed  new  insights  regarding  correlation  between  lattice 
distortion  and  lattice  vibrations  previously  not  possible  owing  to  the  unavailability  of 
experiments.  Now  new  inroads  could  be  made  into  strain  dependent  thermal  engineering  based 
on  the  findings.  We  have  also  performed  first  ever  measurements  of  nanoscale  and  microscale  high 
temperature  creep  in  a  ceramic.  Such  measurements  could  lead  to  significant  advances  in  tunable 
thermal  protection  systems  operating  at  temperatures  ranging  from  very  low  to  ultra-high.  We  have 
used  this  system  for  first  ever  measurements  of  thermal  conductivity  change  in  a  deforming  material 
as  a  function  of  temperature.  Now  experiments  at  elevated  temperatures  up  to  1000  K  are  underway. 
A  request  with  AFoSR  is  pending  to  modify  this  system  to  be  able  to  perform  measurements  at 
temperatures  up  to  2500  K.  Our  experiments  for  the  first  time  show  that  thermal  conductivity  of  even 
a  single  crystal  shows  significant  changes  with  deformation  and  such  deformation  dependent  thermal 
conductivity  is  strongly  correlated  to  nanoscale  stress-strain  curves.  Correspondingly,  thermal 
protection/management  systems  may  perform  in  a  different  way  under  stress  or  environmental 
changes  in  comparison  to  their  original  design.  Ours  is  the  only  experimental  setup  able  to  perform 
such  measurements  at  the  nano  and  micron  scale.  Significant  applications  are  in  the  areas  of  fault 
tolerant  thermal  protection  systems,  topological  meta  materials,  quantum  coherent  sensing  etc..  The 
progress  made  in  this  program  can  also  be  leveraged  for  solving  unique  problems  such  as 
interfacial  thermal  stress/thermal  conduction  problems  at  metal-non/metal  interfaces  (CNT-Au  or 
similar),  innovative  solutions  in  transient  heat  transfer,  and  environment  dependent  heat  transfer 
in  conditions  such  as  nuclear  or  high  temp  environments  or  high  mach  flow  or  corrosive 
environments  etc.  This  framework  can  also  be  leveraged  for  unique  solutions  in  multiphase  or 
small  scale  phase  transformation  heat  transfer  which  has  not  been  possible  till  now  because  of 
limitations  of  classical  molecular  simulations  and  related  experiments. 
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This  research  focused  on  obtaining  an  understanding  of  the  nanoscale  physical  mechanisms  that 
govern  interfacial  interactions  and  heat  transfer  in  materials  at  elevated  temperatures  and  the 
corresponding  thermal  stresses.  The  underlying  strategy  was  to  understand  how  thermal 
conduction  at  nanoscale  gets  affected  by  factors  such  as  the  effect  of  quantum  to  molecular  scale 
temperature  related  phase  change,  effect  of  nanoscale  interfacial  reconstruction,  effect  of  stress 
affected  material  transformation.  Issue  is  to  leverage  this  understanding  to  resolve  thermal  shock 
and  thermal  cycling  affected  failures  in  high  temperature  materials  by  introducing  control  factors 
at  the  nanoscale  which  are  obtained  based  on  understanding  gained  in  this  research.  Material 
systems  of  focus  are  Si  based  ceramics,  particularly  ZrB2-SiC  nanocomposite  interfaces.  For 
verification  purposes  analyses  also  focus  on  Si-Ge  interfaces  and  nanocomposites,  as  a  lot  of 
simulation  data  is  available  on  this  system. 

In  terms  of  modeling  the  required  phenomenon,  classical  analytical  approaches  could  not  be 
used  because  of  small  length  scales  at  material  interfaces.  Classical  molecular  dynamics  or 
Monte  Carlo  methods  based  approach  could  not  be  used  because  classical  interatomic  potential 
used  in  such  approaches  are  not  capable  of  correctly  describing  phase  transformation,  high 
temperature  property  changes,  and  electronic  properties  which  must  be  correctly  described  in 
order  to  properly  understand  such  phenomenon.  There  have  been  significant  advancements  in 
experimental  measurements  of  thermal  conduction.  However,  such  advancements  have  not 
addressed  in-situ  high  temperature  thermal  conduction  as  well  as  thermal  conduction  under 
applied  stress  and/or  strain.  Therefore,  developments  were  made  both  in  terms  of  modeling  and 
experiments  in  this  research  to  model  the  required  material  thermal  behavior.  A  new 
experimental  approach  was  developed  to  perform  in-situ  measurements  of  temperature 
dependent  thermal  conduction  behavior  as  a  function  of  applied  stress/strain.  The  experimental 
developments  were  coupled  with  new  developments  in  non-equilibrium  Green’s  function  based 
ab-initio  calculations  of  nanostructure  thermal  conductivity. 

In  the  following  the  experimental  and  computational  advances  and  the  corresponding  physical 
findings  are  delineated. 

§1  EXPERIMENTAL  DEVELOPMENTS 

Experimental  methods  to  measure  thermal  conductivity  at  microscale  and  nanoscale  structures 
fall  into  two  main  categories  -  electrical  and  optical.  The  difference  between  the  two  types  of 
techniques  is  the  type  of  heating  source  used.  A  range  of  different  shapes  e.g.  thin  films  [1-3],  micro¬ 
cantilevers  [4,  5],  micro-bridges  [6,  7],  and  suspended  membranes  [8-11]  have  been  analyzed  using 
such  techniques.  In  the  electrical  techniques,  metallic  films  deposited  or  machined  on  the  silicon 
structures  are  used  as  a  heating  source  by  applying  a  current  through  the  films  [2].  Temperature  of 
the  silicon  device  at  different  locations  is  measured  by  such  metallic  films  using  electrical-resistance 
thermometry.  In  the  case  of  optical  techniques,  a  focused  laser  source  is  used  as  the  heating  source 
[4],  In  contrast  to  electrical  techniques  optical  ones  are  contact  free.  Therefore,  the  sample 
preparation  is  simpler  than  that  in  the  case  of  other  techniques.  However,  the  difficulty  in  evaluating 
the  absorbance  ratio  of  the  incident  laser  leads  to  uncertainty  in  the  optical  measurements.  Both 
techniques,  however,  can  not  yet  analyze  stress/strain  dependent  thermal  conductivity  as  a  function  in 
a  single  experiment. 


A  new  experimental  method  by  the  Pi’s  group  has  been  established  to  analyze  strain  and 
stress  dependent  thermal  conduction  properties  as  a  function  of  temperature  from  nanoscopic 
length  scale  to  microscopic  length  scale.  In  this  experimental  method  time  resolved  Raman 
Spectroscopy  technique  for  performing  coupled  thermal  and  mechanical  measurements  in  high 
temperature  materials  at  small  scales  has  been  developed.  Another  uniqueness  of  the  setup  is  the 
capability  to  analyze  thermal  and  mechanical  cycling  induced  failure  at  high  temperatures.  At  the 
heart  of  this  setup  is  a  combination  of  a  nanomechanical  stage  with  Raman  spectroscopy,  several 
multispectral  thermometry  stages,  and  a  resistance  thermometry.  The  nanomechanical  stage  also 
incorporates  a  scanning  probe  microscope  to  visualize  the  change  in  material  morphology  under 
straining.  Since  the  chemical  structure  of  the  material  can  only  be  visualized  using  electron 
microscope  observations  (which  is  not  in  the  setup),  quantum  mechanical  simulations  are  used  to 
correlate  scanning  probe  and  Raman  measurements  of  thermal  conductivity  and  mechanical 
deformation.  The  nanomechanical  stage  works  at  a  range  of  temperatures  approaching  750  °C  and 
has  one  nanoscale  and  one  microscale  head  enabling  high  temperature  and  multiscale  measurements. 
During  the  measurements  the  sample  can  be  deformed  to  produce  defects.  In  this  way,  the  impact  of 
defects  on  the  thermal  conductivity  can  be  studied  in  situ.  This  approach  is  inherently  different  from 
the  standard  scenario  where  the  sample  needs  first  to  be  deformed  and  then  studied  for  thermal 
conductivity  afterwards.  In  this  standard  strategy,  the  sample  may  undergo  significant  reconstruction 
when  it  is  transferred  from  the  mechanical  deformation  area  to  the  thermal  measurement.  Thus,  only 
in  situ  measurements  of  deformed  samples  enable  the  real  time  measurement  of  the  impact  of  defect 
formation  effect  on  the  thermal  conductivity.  This  setup  is  uniquely  poised  to  quantitatively 
evaluate  interfacial  stress,  interfacial  thermal  conduction,  and  interfacial  phase  change  as  a 
function  of  temperature  (approaching  1000  K)  in  a  single  in- situ  setting  simultaneously. 


Fig.  1  (a)  Experimental  setup  of  simultaneous  compressing  and  measuring  thermal  conductivity 
as  a  function  of  mechanical  deformation  and  temperature  and  (b)  schematic  of  the  experimental 
arrangement. 
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The  setup  has  been  used  to  study  high  temperature  nanoscale  to  microscale  mechanical 
deformation  in  SiC  based  ceramics  recently.  This  is  first  such  advancement  where  high 
temperature  measurements  upto  500  degree-C  were  performed  for  the  first  time.  Later  on  the 
setup  was  used  to  measure  and  report  thermal  properties  of  a  microscale  Si  Cantilever  as  a 
function  of  applied  stress-strain  and  temperature,  with  such  variations  being  performed  in-situ.  In 
the  following  the  experimental  results  are  described  in  detail. 

§1.1  MEASUREMENT  OF  TEMPERATURE  DEPENDENT  AND  STRESS  DEPENDENT 
THERMAL  PROPERTIES 

The  setup  for  the  measurements  is  shown  in  Fig.  1.  An  AFM  cantilever  (Aspire  conical 
AFM  tip  CT170  from  Nanoscience  Instruments,  Inc.,  AZ)  was  used  as  the  sample.  The  cantilever 
is  made  up  of  highly  doped  single-crystal  silicon.  It  has  rectangular  cross-section  with  the  width 
of  42  /um  and  the  thickness  of  6.5  jum  .  The  length  of  the  cantilever  is  225  jum  .  The  tip  is  located 

at  the  very  end  of  the  cantilever,  where  it  has  a  triangular  shape.  Two  different  sets  of  samples 
made  in  two  different  batches  were  analyzed  (now  referred  to  as  sample  type  1  and  sample  type 
2).  Due  to  the  micrometer  size  of  the  cantilever,  the  heat  flow  through  the  structure  is  low. 
Therefore,  the  precision  of  the  heating  power  measurement  and  temperature  measurement  plays 
an  important  role  in  the  whole  experiment.  In  the  temperature  range  of  our  measurement 
resistance  temperature  detector  (RTD)  is  the  best  choice  for  temperature  measurement.  It  has  a 
tolerance  of  less  than  0.1%,  repeatability  better  than  0.1  °C,  and  response  time  less  than  2s.  The 
standardized  platinum  RTD  devices  have  accurate  resistance-temperature  relationship  after 
calibration,  so  they  are  also  good  candidate  of  precise  heaters.  The  heating  module  consists  of 
two  PT100  RTD  devices  (Minco,  MN,  USA).  It  is  1.7  mm  in  length,  1.3  mm  in  width,  and  0.7 
mm  in  thickness.  In  the  experiment,  a  constant  voltage  was  applied  to  the  heating  circuit.  Both 
the  voltage  and  the  current  (Vi  and  A;  in  Fig.  2,  i=l,2)  were  monitored  as  a  function  of  time.  The 
temperature  at  the  base  end  of  the  cantilever  was  monitored  by  another  RTD  sensor.  Both  RTD 
sensors  were  connected  to  a  RTD  temperature  recorder  (ThermoWorks,  UT,  USA)  with  four- 
wire  bridge  connection  for  the  best  accuracy.  In  the  experiment,  the  silicon  cantilever  was 
mounted  on  a  three-dimensional  linear  stage  with  moving  resolution  of  0.1  pm.  A  resolution  of 
0.1  nm  movement  can  be  achieved  in  the  longitudinal  direction  with  a  piezo-driven  nano¬ 
positioner. 

As  the  length  is  the  dominating  dimension  of  the  cantilever,  it  is  assumed  that  the  heat 
conduction  in  the  longitudinal  direction  dominates  the  overall  heat  flow  from  the  high 
temperature  end  (Th)  to  the  cold  one  (Tc),  if  temperature  gradient  (Th-Tc/Length)  is  maintained 
constant.  The  relation  between  temperature  gradient  and  one  dimensional  steady  state  heat  flow 
is  described  as, 

dT  /IN 

q  =  -kA -  (1) 

dx 

Here,  q  is  the  heat  flow  rate  in  the  cantilever;  k  is  the  thermal  conductivity;  A  is  the  cross- 
sectional  area,  and  dT  /dx  is  the  temperature  gradient  in  the  longitudinal  direction.  The  heat  flow 
rate  q  is  measured  as  an  average  of  the  total  heat  Q  flown  in  a  specific  period.  The  cross-section 
A  is  calculated  from  the  SEM  images  of  cantilever.  dT  (t)/dx  is  derived  from  the  two  RTD 
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sensors.  The  total  energy  flow  Q  is  equal  to  the  difference  between  the  total  input  energy,  Qjn  , 
and  the  energy  dissipated  to  the  surroundings,  Qdis ,  in  this  time  period: 

Q  =  Qin  -  Qdis  (2) 

The  input  energy  Q.n  is  calculated  by  taking  the  integral  of  heating  power  to  time  (Eqn.  3). 


Qin  =  \  "  I'  {t)R(t)dt 


(3) 


The  current  7(f)  flowing  through  the  heater  was  monitored  during  the  experiment.  The 

resistance  of  the  heater  R(t)  is  derived  from  its  resistance-temperature  relationship,  while  its 

temperature  was  read  from  the  RTD  sensor  attached  to  it. 

The  energy  dissipation,  Qdh ,  can  occur  through  lateral  conduction  or  through  convection. 

Because  of  relatively  low  measurement  temperature,  radiation  heat  loss  is  neglected.  Because 
cantilever  length  is  significantly  higher  than  its  cross  sectional  dimensions,  lateral  conduction 
losses  through  cantilever  were  also  neglected.  Dissipated  energy  measurements  by  convection 
and  through  conduction  in  surrounding  support  structures  were  performed  during  each 
measurement  of  thermal  conductivity.  Before  each  experiment,  the  same  voltage  as  intended  for 
thermal  conductivity  measurements  was  applied  to  the  heater  at  zero  load  and  temperature 
fluctuation  was  observed  until  equilibrium  was  attained.  Once  equilibrium  is  attained,  due  to 
steady  state,  the  dissipated  energy  through  convection  is  equal  to  the  input  energy.  It  is  assumed 
that  the  energy  dissipation  rate  remains  constant  during  the  thermal  conductivity  measurements 
that  follow  convective  energy  dissipation  measurements.  This  reasoning  is  based  on  the  fact  that 
during  dissipated  energy  measurements  cantilever  is  kept  2 /um  away  from  the  heater  and  during 

following  measurements  it  is  in  contact  with  the  heater.  The  temperature  change  of  the 
surroundings,  the  heater  and  the  cantilever  during  the  two  different  types  of  measurements  is 
negligible  (<0.5  °C).  Therefore,  its  effect  on  the  energy  dissipation  rate  can  be  safely  assumed  to 
be  negligible. 


Rectangular,  40mN,  9VDC,  400s 


Fig.  2  Temperature  profile  of  the  hot  end  (Th)  and  cold  end  (Tc).  For  the  measurements  shown  a 
constant  load  of  40  mN  was  applied.  The  load  was  held  constant  for  400  s.  The  voltage  of  the 
heating  power  is  9  V  DC. 


11 


Typical  temperature  profiles  of  the  hot  and  cold  ends  are  shown  in  Fig.  2.  Temperature 
profiles  fluctuate  within  a  small  range  (<0.3°C).  This  fluctuation  is  within  1-2  %  of  the  measured 
mean  value,  an  acceptable  error  bound  in  all  experiments  in  the  area  of  reported  experiments. 
The  temperature  drop  from  the  hot  end  to  the  cold  end  is  17  °C.  The  compressive  displacement 
of  the  cantilever  is  derived  from  the  loading  curve  (Fig.  4).  Generally,  the  initial  part  of  the 
loading  curve  has  a  ‘flat  region’,  which  is  caused  by  surface  roughness  of  the  hot  end  heater.  In 
the  analysis,  this  ‘flat  region’  was  removed  by  taking  linear  fitting  of  the  loading  curve. 


Thermal  Conductivity  &  Stress  vs.  Strain 
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Fig.  3.  Thermal  conductivity  as  a  function  of  strain,  stress,  and  temperature  level  with  such 
variations  being  performed  in-situ 


The  calculated  thermal  conductivity  as  a  function  of  applied  strain  and  temperature  is  shown 
in  Fig.  3.  As  shown,  the  thermal  conductivity  shows  an  increase  as  the  compressive  strain  level 
increases  and  decrease  with  increase  in  temperature.  This  trend  of  the  thermal  conductivity 
matches  simulation  results  reported  in  literature,  [13],  for  zero  percent  strain.  The  temperature 
dependence  data  and  stress  dependence  data  is  being  reported  in  this  research  for  the  first  time. 
The  thermal  conductivity  variation  is  as  large  as  20  %  to  40  %  depending  upon  the  applied  strain. 
The  corresponding  stress  percentage  change  is  almost  the  same  as  the  thermal  conductivity 
percentage  change.  This  finding  indicates  that  both  thermal  conductivity  and  stress  in  a  sample 
can  vary  in  the  similar  manner  as  a  function  of  applied  strain.  This  is  particularly  important  for 
small  scale  devices  which  operate  at  high  number  densities  and  at  high  thermal  stress  levels. 
Thermal  management  analyses  of  such  devices  based  on  0%  strain  thermal  conductivity  values 
may  give  incorrect  prediction  of  thermal  stress. 

The  thermal  conductivity  measurements  at  elevated  temperatures  approaching  1000  K  are  in 
progress.  Such  setup  combines  the  results  just  described  with  author’s  work  on  high  temperature 
nanoscopic  and  microscopic  mechanical  property  measurements  which  were  partially  supported 
by  this  project  funding.  The  next  section  describes  and  discusses  findings  in  those  measurements. 


§1.2  MEASUREMENT  OF  HIGH  TEMPERATURE  PROPERTIES 

Nanoindentation  creep  has  been  observed  in  a  wide  range  of  materials  including  glasses, 
ceramics,  and  metallic  materials,  [12-14],  The  indentation  creep  rate  for  high  melting  point 
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materials  can  be  of  the  order  of  5  x  10~5  s'1.  This  is  a  large  value  when  compared  to  the 
corresponding  bulk  strain  rate.  The  models  for  indentation  creep  are  the  same  as  those  for  bulk 
creep,  with  the  equivalent  stress  exponent,  n,  used  as  an  indicator  of  the  creep  mechanism,  [12],  It 
is  believed  that  when  the  value  of  n  is  1,  creep  is  controlled  by  vacancy  diffusion  as  deformation 
mechanism,  [15,  16];  when  the  n  value  is  2,  the  creep  mechanism  the  controlled  by  grain 
boundary  sliding,  [17];  when  n  is  3,  diffusion-controlled  dislocation  motion  dominates  as 
deformation  mechanism,  [18,  19];  and  when  n  is  5,  it  is  dislocation  climb-controlled  creep 
mechanism,  [20].  During  indentation  creep  tests  on  certain  metals,  alloys,  and  ceramics  at  room 
temperature,  high  n  values  up  to  hundreds  have  been  observed,  [14,  21-26].  The  mechanism 
behind  such  high  stress  exponent  values  has  been  attributed  to  volumetric  densification  and 
dislocation  pile  up.  Materials  systems  analyzed  include  metallic  thin  films,  e.g.  Cu,  Al,  Ni,  [27- 
29],  thin  films  made  of  metallic  alloys,  e.g.  Ni3Al,  [30],  and  ceramics  e.g.  Zr02,  [31],  and  SiC- 
based,  [32-34].  During  such  studies,  the  size  of  the  indenter  (owing  to  indentation  size  effect- 
ISE)  and  the  grain  size  of  the  polycrystalline  materials  have  been  shown  to  strongly  affect  the 
indentation  creep  response  of  materials.  Besides  the  indentation  size  effect  (ISE)  and  the  grain 
size  effect,  the  creep  deformation  mechanism  is  also  shown  to  be  affected  by  the  test 
temperatures  and  the  change  in  microstructure  of  the  material  as  a  result  of  temperature  change. 
It  has  also  been  shown  that  fiber  reinforcement,  [35],  or  particle  reinforcement,  [36],  can 
constrain  creep  rate.  This  work  for  the  first  time  presented  nanoindentation  and  microindentation 
creep  analyses  on  polymer  derived  Si-C-O  ceramic  coatings  at  temperatures  ranging  from  room 
temperature  to  500  degree-C. 


(a)  (b) 

Fig.  4:  (a)  Nanoindentation  elastic  moduli  and  (b)  microindentation  elastic  moduli  of  Si-C-O 
coatings  as  a  function  of  temperature  and  maximum  indentation  load 


Fig.  4  shows  the  elastic  moduli  values  as  a  function  of  testing  temperature  and  peak 
indentation  load.  As  shown,  the  elastic  moduli  increase  with  increase  in  temperature  and  with 
increase  in  peak  load  at  the  nanoscale.  At  the  microscale  the  elastic  moduli  increase  as  a  function 
of  temperature.  However,  dependence  on  the  peak  load  at  the  microscale  is  not  clear.  Increase  as 
a  function  of  peak  load  signifies  the  ISE.  It  has  been  shown  that  Young’s  modulus  of  SiC>2  will 
increase  with  temperature  due  to  compaction  and  distortion  of  Si04  tetrahedra  in  the  temperature 
range  of  experiments.  The  underlying  SiCO  matrix  retains  this  attribute.  An  examination  of  Fig. 
4  (a)  reveals  the  ISE  is  stronger  at  higher  temperatures.  The  ISE  is  absent  in  the  case  of 
microscale  data.  The  reason  behind  the  observation  of  ISE  at  the  nanoscale  is  not  clear.  For  each 
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data  point  in  Fig.  4  (a),  a  minimum  of  16  tests  were  performed.  Therefore  a  possibility  that  a 
single  microstructural  aberration  is  contributing  to  the  ISE  is  ruled  out.  At  the  nanoscale  the 
nanoindenter  radius  is  approximately  20  nm  and  projected  area  is  approximately  0.5  pm2. 
Accordingly,  the  indenter  samples  either  the  SiC  particles  or  small  spaces  of  SiCO  matrix 
between  adjacent  TiSi2  particles.  The  ISE,  therefore,  can  be  attributed  to  increasing  nanoscale 
densification  with  increasing  indentation  depth.  At  the  microscale  collective  deformation  of  a 
sample  volume  containing  many  SiC  particles  takes  place.  Due  to  such  collective  motion  of 
particles  under  indenter  tip,  the  effect  of  volumetric  densification  under  indenter  tip  is  minimzed 
leading  to  the  absence  of  ISE.  Increase  in  temperature  contributes  to  increased  binding  of  TiSi2- 
SiCO  interfaces  invariably  leading  to  increase  in  the  elastic  moduli  values  as  a  function  of 
temperature  at  both  scales. 
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Fig.  5:  (a)  Nanoindentation  hardness  and  (b)  microindentation  hardness  of  Si-C-O  coatings  as  a 
function  of  temperature  and  maximum  indentation  load 


Figure  5  plots  the  hardness  data  corresponding  to  the  moduli  plots  in  Fig.  5.  As  shown,  at  the 
nanoscale  the  trend  for  the  elastic  moduli  is  repeated  for  the  hardness  values.  A  clear  trend  on 
increase  in  Meyer’s  exponent  with  increase  in  temperature  emerges  from  nanoindentation  data. 
The  microindentation  hardness  data  trend  is  opposite  to  the  trend  shown  in  the  case  of 
microindentation  elastic  moduli.  Hardness  decreases  as  the  measurement  temperature  increases. 
In  addition,  hardness  reduces  with  increase  in  the  peak  indentation  load,  signifying  strain 
softening  of  the  material.  Reduction  of  hardness  with  increase  in  temperature  is  attributable  to 
stronger  TiSi2  particles  pressing  in  relatively  softer  SiCO\matrix  getting  progressively  softer  with 
increase  in  temperature  at  the  microscale.  At  the  nanoscale  not  enough  contact  area  is  available 
to  have  such  an  effect.  Data  on  correlation  between  temperature  dependent  hardness  and  young’s 
moduli  values  for  ceramics  is  not  available  in  the  literature.  A  similar  trend  (increase  in  moduli 
and  decrease  in  hardness  as  a  function  of  temperature)  has  been  observed  for  Chromium  steels 
by  Medved  and  Bryukhanov  [37], 

It  has  been  shown  that  in  Si-C-0  PDCs  in  amorphous  phase  the  deformation  during 
indentation  occurs  due  to  volumetric  densification  events.  It  has  been  shown  that  volumetric 
densification  promotes  strain  hardening  in  vitreous  silica  under  both  indentation  and  diamond- 
anvil  compression  experiments,  [34],  The  microindentation  experiments  on  the  present  materials 
in  the  amorphous  state  reveal  a  reduction  of  indentation  hardness  under  increasing  applied  loads, 
indicating  densification-assisted  strain  softening  of  these  materials,  [34].  A  combination  of 
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nanoscale  and  microscale  results  points  that  individually  neither  the  SiCO  matrix  nor  TiSi2 
particles  show  strain  softening.  However,  mechanistic  effects  related  to  their  combined  presence 
at  the  microscale  results  in  strain  softening  behavior.  Cheng  and  Cheng,  [38],  have  pointed  out 
that  indentation  hardness  is  not  a  function  of  indentation  depth  for  materials  that  truly  follow 
power-law  strain  hardening.  This  is  in  contrast  to  the  nanoindentation  and  microindentation 
results  from  the  present  materials  that  show  an  increase  in  hardness  with  increasing  load  during 
nanoindentation  and  reduction  in  hardness  with  increasing  load  during  microindentation.  Apart 
from  the  influence  of  strain  hardening  or  softening,  material  pile-up  around  the  indent  and 
indentation  creep  can  also  contribute  to  indentation  hardness  behavior  observed  in  the  present 
work. 
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Fig.  6:  Indentation  depth  (original  and  thermal  drift  corrected)  as  a  function  of  time  at  different 
temperatures  in  the  case  of  (a)  nanoindentation  test  at  the  peak  load  of  300  mN,  (b) 
nanoindentation  test  at  the  peak  load  of  500  mN,  (c)  microindentation  test  at  the  peak  load  of  750 
mN,  and  (d)  microindentation  test  at  the  peak  load  of  1000  mN 


Figure  6  plots  creep  depth  as  a  function  of  dwell  period  at  three  different  temperatures  at  two 
different  nanoindentation  and  microindentation  peak  loads.  Creep  data  is  normalized  by 
subtracting  the  initial  depth  of  each  creep  test.  For  comparison,  the  plots  also  show  the  thermal 
drift  uncorrected  as  well  as  corrected  data.  As  expected,  the  effect  of  thermal  drift  correction  is 
the  highest  at  the  highest  temperature  and  at  the  nanoscale.  The  creep  depth  vs.  time  plots 
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reached  steady  state  within  first  100  s  of  plotting  in  all  cases.  Effect  of  temperature  is  to  increase 
the  creep  rate.  The  temperature  effect  is  more  pronounced  at  the  micro  scale  in  comparison  to  the 
nanoscale.  The  trend  is  particularly  strong  when  transitioning  from  250  °C  to  500  °C. 


Nano-Creep  of  Si-C-O,  Thermal  Drift  Corrected 
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Fig.  7:  (a)  stress  exponent  as  a  function  of  load  during  nanoindentation  creep  measurements,  (b) 
stress  exponent  as  a  function  of  load  during  microindentation  creep  measurements,  (c)  creep 
strain  rate  as  a  function  of  load  during  nanoindentation  creep  measurements  and  (d)  creep  strain 
rate  as  a  function  of  load  during  microindentation  creep  measurements 


Figure  7  plots  stress  exponent  and  creep  strain  rate  as  a  function  of  temperature  and  peak 
indentation  load  at  both  length  scales.  Stress  exponent  at  the  nanoscale  lies  in  the  range  of  4  to  5 
indicating  a  dislocation  climb  related  creep  deformation  mechanism  occurring  through  bulk  or 
pipeline  diffusion  of  dislocations.  The  exponent  reduces  with  increase  in  temperature  indicating 
a  transition  of  mechanism  from  dislocation  climb  to  diffusion.  The  densely  dislocated  structure 
with  high  dislocation  density  under  indenter  may  result  in  dense  dislocation  channels  directing 
atoms  from  dislocated  cores  to  free  surfaces  resulting  in  such  a  mechanism  occurring  at  low 
homologous  temperatures.  At  the  microscale  the  stress  exponent  is  considerably  higher 
indicating  that  primary  mechanism  of  deformation  is  volumetric  densification.  Such  a  transition 
suggests  a  rapid  change  of  mechanism  from  linear  diffusion  flow  to  power  law  mechanism  (e.g. 
climb)  and  eventually  to  rate  insensitive  plastic  flow  (dislocation  glide)  as  the  indentation  size 
moves  towards  microscale  regime.  The  disparity  in  stress  exponent  values  at  the  two  length 
scales  is  not  reflected  in  the  creep  strain  rate  data.  As  shown,  the  creep  strain  rates  at  both  length 
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scales  lie  in  the  same  range  with  microscale  data  being  on  the  lower  side.  Creep  strain  rate 
generally  increases  as  a  function  of  increase  in  temperature  as  well  as  with  increase  in  peak 
indentation  load.  The  creep  strain  rate  at  the  microscale  being  lower  than  that  at  the  nanoscale 
can  be  attributed  to  the  averaging  effect  of  measurements.  As  discussed  earlier,  in  the  case  of 
microscale  measurements,  the  indenter  covers  an  area  that  is  approximately  10  times  larger  than 
the  nanoindenter  area.  Since  hardness  reduces  and  creep  strain  rate  increases  with  increase  in 
temperature  at  the  microscale,  the  strain  rate  sensitivity  index  (constant  k  in  relation  a=b  s  k) 
indicates  strain  softening  at  the  microscale.  The  same  analyses  reveal  strain  hardening  at  the 
nanoscale. 
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Fig.  8:  Thermal  activation  volume  as  a  function  of  peak  indentation  load  and  measurement 

temperature  at  (a)  nanoscale  and  (b)  microscale. 


In  general,  in  the  case  of  PDCs  both  dislocation  dominated  and  volumetric  densifying  type  of 
deformation  mechanisms  with  varying  degree  are  present  during  indentation.  In  the  case  of 
quartz,  the  size  effect  on  the  stress  exponent  can  be  explained  by  a  reduction  of  the  localized 
shear  volume  as  the  indent  size  decreases.  However,  in  the  present  work  the  size  effect  is  related 
to  shift  from  dislocation  climb  and  diffusion  related  mechanism  at  the  nanoscale  to  volumetric 
densification  at  the  microscale  (Fig.  8).  The  extent  of  a  particular  type  of  mechanism  determining 
deformation  is  governed  by  the  availability  of  excess  free  volume  which  is  related  to  thermal 
activation  volume,  which  is  also  the  average  change  in  volume  of  the  flow  unit.  Low-temperature 
creep  deformation  strain  rate,  s  ,  of  materials  can  be  expressed  by  Arrhenius-type  flow  function, 
[39], 
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In  the  above  Eqn.,  e  is  reference  strain  rate  that  varies  with  the  type  of  material,  stress  level, 


and  microstructure,  ag  is  the  activation  energy  for  creep  or  other  rate-dependent  process,  kB  is 
Boltzmann’s  constant  and  T  is  temperature  in  Kelvin.  Parameter  e  is  proportional  to  the 


concentration  of  elementary  defects  which  cause  plastic  strain,  [40],  The  sensitivity  of  strain  rate 
to  stress  is  mainly  determined  by  the  a  g  term.  The  thermal  activation  volume  (v  * )  is  expressed 
by  the  partial  derivative  of  ag  with  respect  to  stress  a  ,  [41].  In  the  present  work,  v  *  is 
approximated  as 
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The  derivatives  are  taken  at  fixed  strain  at  each  temperature  to  constrain  the  influence  of 
structure  change  due  to  work  hardening  and  softening.  Figure  10  plots  the  thermal  activation 
volume  calculated  in  this  case  as  a  function  of  temperature,  peak  indentation  load,  and  length 
scale.  As  shown,  thermal  activation  volume  increases  with  increase  in  the  length  scale  of 
measurements.  Since,  the  thermal  activation  volume  is  directly  correlated  to  the  deformation 
energy,  the  trend  implies  a  significantly  high  deformation  work  required  at  the  microscale  in 
comparison  to  the  nanoscale.  An  almost  linear  increase  in  activation  volume  with  increase  in 
temperature  implies  that  deformation  mechanism  is  primarily  controlled  by  dislocations.  At  the 
nanoscale  the  thermal  activation  volume  is  10  times  less  than  that  at  the  microscale  indicating 
that  the  dislocation  pile  up  is  restricted  as  well  as  that  dislocation  climb  and  diffusion  is  the  main 
deformation  mechanism.  At  microscale  significantly  higher  thermal  activation  volume  indicates 
availability  of  excess  free  volume  in  the  material  structure  that  can  be  eliminated  by  densification. 
A  combination  of  earlier  findings  with  this  result  indicates  that  with  increase  in  length  scale  the 
deformation  mechanism  switches  from  dislocation  diffusion  to  dislocation  pileup.  With  increase 
in  length  scale,  dislocation  motion  does  not  face  significant  obstacles  or  phase  transformation 
related  factors.  Such  factors,  if  present,  will  change  the  almost  linear  trend  observed  in  Fig.  10. 
The  observed  trend  of  the  data  on  peak  indentation  load  as  well  as  stress  exponent  values  at  room 
temperature  are  in  close  agreement  with  the  work  of  Li  and  Ngan  [14].  However,  stress 
exponents  are  lower  than  those  reported  by  Li  and  Ngan  for  fused  silica,  indicating  a  difference 
in  deformation  mechanism.  The  deformation  mechanism  is  intricately  linked  with  the 
microstructure  of  this  material.  The  material  has  TiSi2  particles  with  higher  melting  point 
embedded  in  SiCO  matrix  with  lower  melting  point.  With  increase  in  temperature,  softening  of 
the  matrix,  particle  sliding,  and  particle  rearrangement  contribute  to  deformation  mechanism. 
The  effect  of  these  factors  is  more  pronounced  at  microscale  due  to  higher  surface  area  sampled. 

§2  COMPUTATIONAL  DEVELOPMENTS 

Computational  advances  focused  on  establishing  an  ab-initio  approach  to  predict  chemical 
composition  and  microstructure  dependent  thermal  conductivity  of  examined  materials  with  an 
account  of  nanoscale  features.  As  pointed  out  earlier,  classical  molecular  simulations  are  limited 
in  their  applicability  to  predict  thermal  properties  of  materials  at  high  temperatures  under  rapid 
transient  heating  due  to  transient  microstructure  and  chemical  evolution  occurring  during  such 
imposed  conditions.  Classical  molecular  simulations  are  also  limited  only  to  analyzing  thermal 
transport  in  phonon  dominated  systems.  A  widely  studied  semiconductor  such  as  Si  becomes 
conductor  at  high  temperatures  (ban-gap  reduces  as  a  function  of  temperature  increase)  severely 
limiting  classical  molecular  simulation  approach’s  applicability  to  analyze  small  scale  heat 
transfer  even  in  such  a  material.  In  this  research,  a  new  quantum  mechanical  method  based  on 
extension  of  non-equilibrium  Green’s  function  method  is  established  that  could  be  used  to  study 
fundamental  issues  affecting  phase  change  heat  transfer  at  quantum  and  molecular  scales. 
Particular  application  relevant  to  project  is  being  able  to  resolve  issues  regarding  grain  boundary 
and  interface  composition  and  how  phase  change  and  stresses  at  high  temperature  affect  thermal 
properties  that  lead  to  thermal  shock  failure.  This  method  is  used  to  study  interfacial  heat  transfer 
in  ZrB2-SiC  material  system  of  interest  to  AFoSR.  Since  ab-initio  approach  is  limited  to  small 
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system  sizes  because  of  significant  computational  demand,  such  simulations  are  being  coupled  to 
approximate  classical  molecular  simulations  to  scale  up  the  simulation  predictions  to  realistic 
microstructure  size.  Accordingly,  the  development  of  ab-initio  approach  was  performed  together 
with  analyses  of  interfacial  heat  transfer  as  a  function  of  molecular  scale  morphology.  Analyses 
have  proven  for  the  first  time  that  materials  with  biomimetic  phase  morphology  have  thermal 
conductivity  values  independent  of  strain.  This  finding  has  strong  implication  for  developing 
materials  with  thermal  properties  independent  of  applied  stress.  We  have  proven  for  the  first  time 
that  that  tensile  straining  and  heat  flow  direction  can  be  used  to  develop  a  thermal  diode  material 
from  superlattice  construction.  In  addition,  our  group  was  the  first  to  show  that  nanostructures 
with  tunable  thermal  properties  could  be  developed  based  on  strain  engineering.  For  this  reason, 
the  relevant  publication  was  featured  in  the  Virtual  J.  of  Nanoscale  Science  and  Technology  (A 
collection  of  significant  advances  in  nanotechnology).  In  the  following  both,  the  ab-initio 
calculations  and  the  classical  molecular  simulations  are  briefly  presented. 

§2.1  NON-EQUILIBRIUM  GREEN’S  FUNCTION  BASED  FORMALISM  TO  STUDY 
INTERFACIAL  HEAT  TRANSFER  AT  HIGH  TEMPERATURES  AND/OR  UNDER 
TRANSIENT  HEAT  LOAD  CONDITIONS 

Though  thermal  conductance  of  a  material  depends  both  on  phonons  and  electrons, 
significant  percentage  of  past  research  has  primarily  addressed  the  phononic  contributions.  With 
change  in  temperature  the  ratio  of  phonon  thermal  conductance  to  electron  thermal  conductance 
changes.  It  has  been  argued  in  ceramics  community  that  under  transient  heat  loads  such  as  those 
during  thermal  shock  loading,  electron  thermal  conductance  can  play  an  important  role  in  a 
material’s  response.  Therefore,  understanding  phonon  and  electron  contributions  to  a  material’s 
thermal  response,  especially  at  high  temperatures,  is  an  important  issue. 

Phononic  contribution  to  the  thermal  conductance  has  been  traditionally  analyzed  using 
modeling  approaches  based  on  the  Boltzmann  Transport  Equation  (BTE).  In  the  instances  where 
an  interface  is  being  analyzed,  variants  of  either  acoustic  mismatch  model  (AMM)  or  of  diffusive 
mismatch  model  (DMM)  are  employed  for  the  description  of  conductive  heat  transport.  The 
AMM  makes  the  assumption  of  thermal  equilibrium  achieved  immediately  on  either  side  of  a 
single  abrupt  junction.  Calculations  of  the  AMM  only  consider  the  materials  on  either  side  of  the 
interfacial  bond  and  not  the  interfacial  bond  itself.  The  bond  at  the  atomic  interface,  however, 
has  been  experimentally  shown  to  drive  interfacial  transport  and  scattering,  [42],  DMM  approach 
too  has  been  shown  to  be  limited  in  effectively  capturing  the  conductive  heat  transport  across  the 
interfacial  zone.  There  is  no  clear  recipe  on  how  to  combine  individual  impedances  for  an 
extended  channel  region  driven  into  non-equilibrium  by  two  separate  thermal  reservoirs,  as  is 
commonly  invoked  for  electron  flow.  Both  methods  ignore  evolution  of  nanostructure  at 
interfaces  as  a  function  of  temperature  or  applied  strain  due  mainly  to  limitation  of  classical 
interatomic  potentials  in  ability  to  capture  complex  phase  changes  that  may  be  occurring  at 
interfaces,  especially  at  high  temperatures.  Both  methods  also  cannot  account  for  electron 
thermal  conductance.  Electron  flow  has  usually  been  formulated  as  a  two-junction  problem, 
making  a  clear  conceptual  separation  between  “contacts”  that  are  held  at  separate  thermal 
equilibrium  by  a  battery,  aided  by  inelastic  thermalizing  scattering  events  in  them,  and  a 
“channel”  or  “device”  region  where  the  actual  transport  dynamics  occurs.  The  viewpoint  leads 
naturally  to  the  non  equilibrium  Green’s  function  (NEGF)  technique.  The  NEGF  formalism 
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calculates  wave  (electron  or  phonon)  transport  across  the  channel,  and  therefore  takes  into 
account  masses  and  forces  bond  strengths  on  either  side  of  the  channel  atom. 

The  phononic  thermal  conductance  calculation  in  the  present  work  is  based  on  the  phonon 
density  of  states  and  the  frequency  calculations.  Electronic  part  of  thermal  conductance  is 
calculated  using  a  recent  modification  to  the  quantum  mechanical  NEGF  method,  [43],  The 
NEGF  approach  had  its  origins  in  the  development  of  the  theory  of  quantum  electrodynamics  in 
the  late  1940's  and  early  1950's,  and  inherits  this  field's  emphasis  on  perturbation  expansions 
described  in  diagrammatic  form.  The  NEGF  theory  was  formulated  by  Kadanoff  et  al  [44]  and 
Keldysh  [45],  It  was  later  elaborated  by  Langreth  [46]  as  described  in  the  text  by  Mahan  [47], 
The  problems  and  promise  of  applying  this  approach  to  electron  devices  has  been  discussed  by 
Huag  et  al.  [48].  NEGF  approach  has  been  extensively  used  to  study  electronic  transport  in 
nanostructures  in  the  electronics  domain,  [49]  and  [48].  The  analyses  regarding  the  role  of 
electrons  and  how  their  behavior  changes  temperature  dependent  properties  have  been  limited  to 
a  certain  class  of  material  systems  due  to  extremely  high  computational  demands  of  NEGF 
method.  A  change  was  brought  by  Wang  et  al.  [50]  who  applied  conventional  electronics- 
focused  approach  to  calculate  electronic  thermal  conductance  of  single  chain  semiconductor 
atoms.  They  proposed  significant  modifications  in  the  Feynman  Diagrams  to  simplify  it  to  be 
applicable  to  systems  with  high  number  of  electrons.  An  alternate  approach  is  based  on  atomistic 
Green’s  function  calculation  by  Zhang  et  al.  [51],  However,  in  this  approach  one  uses  classical 
interatomic  potentials  and  therefore  cannot  account  for  either  complex  atomistic  nanostructure 
evolution  at  an  interface  under  complex  tri-axial  stress  states,  especially  at  high  temperatures  or 
electron  mediated  as  well  as  electron-phonon  coupling  mediated  conductive  thermal  transport. 
Advancement  is  needed  to  address  these  issues  so  that  conductive  heat  transport  in  realistic 
conditions  in  a  material  interface  can  be  examined.  Wang  et  al.  [52],  generalized  the  NEGF 
approach  by  replacing  fermions  with  bosons  as  basic  entities.  However,  such  analyses  have  been 
limited  to  1-D  atomic  chains  due  to  the  computational  demand  and  extreme  complexity  of  the 
method.  The  applicability  of  modified  NEGF  to  boson  based  analysis  was  further  utilized  by 
various  groups  such  as  pointed  out  in  refs  [53],  [54],  and  [55].  But  this  too  has  been  limited  to  1- 
D  atomic  chains  due  to  the  computational  demand  and  extreme  complexity  of  the  method. 

In  this  work,  NEGF  is  modified  to  make  it  applicable  to  bigger  atomistic  systems  wherever 
one  could  to  strike  a  balance  between  accuracy  of  the  calculations  and  reduction  in  the 
complexity  and  computational  demand,  [43].  This  approach  provides  us  with  the  capability  to 
separate  out  the  phononic  and  electronic  contribution  to  the  overall  thermal  conductance.  One 
significant  aspect  of  the  current  work  is  that  the  whole  framework  is  based  on  first  principles 
calculations  and  does  not  use  a  classical  interatomic  potential  like  those  used  in  atomistic 
Green’s  function  formulations.  The  use  of  first  principles  framework  enables  applicability  of 
current  work  in  complex  nanoscale  material  situations  such  as  phase  change,  high  temperature 
behavior  etc.  The  need  for  any  interatomic  potential  functions  has  been  removed  by  calculating 
the  parameters  based  on  purely  electronic  and  phononic  eigen  function  analysis.  Such  use  enable 
addressing  complex  atomic  structures  that  may  also  involve  phase  transformation. 

The  NEGF  approach  for  phonon  transport  is  implemented  in  two  steps:  first  the  phonon 
propagation  equations  for  the  contacts  are  partitioned  out  to  create  an  equivalent  open-boundary 
version  of  Newton’s  law  for  lattice  vibrations  in  the  channel.  The  contact  states  enter  their 
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boundary  condition  through  carefully  computed  self-energy  matrices.  Next,  a  thermal 
equilibrium  condition  is  imposed  on  the  contact  state  bilinear  variables  in  order  to  compute  the 
NEGF  thermal  current.  The  NEGF  method  has  as  its  main  ingredient  in  the  Green’s  function, 
which  is  a  function  of  space-time  coordinates.  From  knowledge  of  this  function  one  can  calculate 
time-dependent  expectation  values  such  as  currents  and  densities,  electron  addition  and  removal 
energies  and  the  total  energy  of  the  system.  NEGF  can  be  applied  to  both  extended  and  finite 
systems.  To  describe  many-electron  systems  in  general  time-dependent  external  fields  NEGF 
approach  here  restricts  itself  to  a  non-relativistic  description  based  on  the  time-dependent 
Schrodinger  equation  (TDSE).  We  are  utilizing  this  simplified  NEGF  approach  based  on  first 
principles.  The  thermal  conductance  of  the  material  system  can  be  calculated  by  the  expression 
provided  by  Wang  et  al.  [52]  also  known  as  Taylor  expansion  of  the  Landauer  Formula  for 
thermal  conductance, 
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Here,  g  '  is  the  retarded  Green's  Function  while  G  “  is  the  advanced  Green's  Function,  r  L  and 
r  n  represent  the  linear  and  non-linear  part  of  the  frequency  broadening  terms,  and  S  represents 
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the  interaction  term.  - is  the  probability  distribution  function  given  by  Bose-Einstein 
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)  in  the  case  of  phononic  contribution  and  Fermi-Dirac 
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)  in  the  case  of  electronic  contribution.  Note  that  although 


the  ingredients  in  this  equation  are  calculated  using  open  boundary  Newton’s  law,  the 
conductance  equation  itself  is  quantum  mechanical,  as  it  apportions  a  fixed  phase  space  for  each 
thermal  mode. 

For  analyzing  interfacial  heat  transfer  SiC[lll]-ZrB2[0001]  was  chosen  as  the  model  system. 
This  interface  was  studied  along  with  single  crystalline  SiC  and  ZrB2  systems.  The  interface  is 
shown  in  Fig.  9.  Since  substantial  temperature  and  strain  related  simulations  along  with 
significant  post  processing  calculations  need  to  be  performed,  only  one  interface  is  analyzed  due 
to  computational  time  constraints  (simulations  reported  took  1  and  a  half  years).  All  supercells 
have  periodic  boundary  conditions  imposed  in  all  directions.  Interfacial  supercell  contains 
sufficient  vacuum  space  above  and  below  the  atomic  slabs  to  minimize  coupling  between  their 

O 

free  surfaces.  Because  of  the  large  lattice  mismatch  between  ZrB2  (<7zrB2  =3.169  A)  and  SiC 


21 


(< dsicim ]  =3.082  A)  we  follow  approach  of  Trivedi  et  al.  [56]  and  adopt  cell  dimensions  which 
impose  a  tensile  strain  on  SiC  to  accommodate  the  stiffer  and  considerably  thicker  (bulk-like) 
ZrB2  template.  We  further  assume  that  the  large  strain  differential  associated  with  this  coherent 
interface  approximation  does  not  alter  the  energetics  of  the  interfacial  bonding. 


Fig.  9  The  interface  examined  in  the  present  research.  The  single  crystal  samples  have  the  same 
overall  size  and  orientation  with  respect  to  stretching  direction  depicted  using  arrows. 

Structure  calculations  of  the  SiC  [lll]-ZrB2  [0001]  heterojunction  were  carried  out  using  the 
NWChem  simulation  package,  [57].  This  program  employs  pseudopotentials  derived  from  the 
plane- wave  formalism  and  the  Perdew91  form  of  the  GGA  functional.  We  used  standard  norm- 
conserving  Trouiller-Martin’s  and  Hamann  pseudopotentials.  A  plane  wave  cutoff  energy  of  400 
eV  was  used  in  conjunction  with  a  real  space  grid  of  16x16x16  to  achieve  a  force  accuracy  of 

O 

0.01  eV/A.  This  grid  size  was  chosen  after  studying  different  grid  sizes  for  force  convergence. 
The  atomic  coordinates  were  fully  relaxed  to  their  zero  force  positions,  yielding  optimized 
atomic  structures. 


Fig.  10  (a)  Tensile  peak  strength  of  single  crystal  ZrB2,  single  crystal  SiC,  and  ZrB2-SiC 
interface  as  a  function  of  temperature  and  (b)  average  interatomic  distance  in  single  crystal  SiC, 
ZrB2,  and  SiC-ZrB2  interface  as  a  function  of  temperature.  This  plot  also  shows  experimental 
measurement  of  average  SiC  interatomic  distance  reported  by  Fahrenholtz  et  al.  (2007)  as  a 
function  of  temperature.  This  is  the  only  such  measurement  available. 

Figure  10-(a)  shows  peak  strength  as  a  function  of  temperature  for  the  examined  SiC  and 
ZrB2  single  crystals  and  the  SiC-ZrB2  interface.  The  strength  is  calculated  as  the  peak  value  of 
stress  as  all  the  examined  supercells  are  stretched.  The  interface  is  stretched  in  direction  shown 
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using  arrows  in  Fig.  9  earlier.  As  shown  the  strength  of  all  samples  reduces  with  increase  in 
temperature.  Strength  of  interface  is  higher  than  that  of  ZrB2  and  SiC  single  crystals,  as  expected. 
It  is  important  to  note  here  that  in  the  examined  structures  structural  changes  were  observed  as  a 
function  of  strain  that  include  bond  rotation  and  possibly  defect  formation.  Due  to  complex 
crystal  structure  it  is  not  possible  to  explicitly  characterize  the  type  of  defects  formed  (screw  vs. 
edge  dislocation,  the  extent  of  individual  bond  rotations  due  to  a  large  number  of  bond  angle 
permutations  involved  etc.).  One  way  to  characterize  strength  loss  in  the  samples  as  a  function  of 
temperature  is  to  compare  average  bond  length  of  the  samples  as  a  function  of  temperature. 
Figure  10-(b)  shows  average  bond  length  in  the  systems  examined  as  a  function  of  temperature. 
As  shown  in  Fig.  10-(b)  average  bond  length  increase  in  all  three  samples  can  be  correlated  to 
the  loss  of  strength  as  a  function  of  temperature.  Change  in  average  bond  length  as  a  function  of 
temperature  also  showed  a  similar  trend  (not  shown).  However,  the  change  is  a  very  small 
fraction  (2-5%)  of  the  average  bond  length  values  shown  indicating  that  straining  has  significant 
bond  rotations  involved.  Overall,  however,  it  is  interesting  to  note  that  average  bond  length  value 
correlated  proportionally  to  the  loss  of  strength  as  a  function  of  temperature.  The  average  bond 
length  change  data  as  a  function  of  temperature  change  is  not  available  in  literature.  However, 
Fahrenholtz  et  al.  [58]  recently  have  examined  average  bond  length  as  a  function  of  temperature 
in  SiC  phase  of  the  ZrB2-SiC  composites  prepared  in  their  lab.  The  trend  for  SiC  matches  with 
that  reported  by  Fahrenholtz  et  al.  [58],  It  is  important  to  note  that  SiC  data,  here,  is  for  the  single 
crystal  and  Ref  [58]  data  is  for  SiC  phase  in  a  composite. 


SiC-ZrB2  interface  as  a  function  of  strain  and  temperature  and  (b)  percentage  change  with 
respect  to  500  K  value  in  phonon  thermal  diffusivity  in  single  crystalline  SiC,  single  crystalline 
ZrB2,  and  SiC-ZrB2  interface  as  a  function  of  strain  and  temperature 

Fig.  11 -(a)  shows  phononic  contribution  to  thermal  conductance  as  a  function  of  temperature 
and  at  0%  and  10%  strain  values.  As  shown,  phonon  thermal  conductance  of  SiC  reduces  as  a 
function  of  temperature  at  both  0%  and  10%  strain  values.  However,  in  the  case  of  ZrB2  the 
phononic  contribution  to  thermal  conductance  increases  as  function  of  temperature.  This  can 
possibly  be  attributed  to  combined  sp-sp2-sp 3  hybridized  state  of  bonds  in  ZrB2  (besides  highly 
polarized  d  orbitals  which  also  lead  to  spd  hybrid  bond  configuration)  instead  of  the  sp3- 
hybridized  bonds  in  the  case  of  SiC.  The  slightly  metallic  character  of  bonds  in  ZrB2  leads  to  free 
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electrons  and  slightly  ionic  cores  participating  in  thermal  conduction  as  against  in  the  case  of  SiC 
where  bond  character  is  strongly  covalent  with  no  electrons  available  for  thermal  conduction. 
With  increase  in  temperature  the  thermal  energy  in  ZrB2  gets  distributed  between  ionic  cores  and 
free  electrons.  Increased  mobility  of  free  electrons  at  increasing  temperature  also  contributes  to 
ionic  core  mobility  with  the  remaining  thermal  energy  getting  deposited  in  the  acoustic  phonon 
mods  leading  to  increase  in  phonon  thermal  conductance.  The  thermal  diffusivity  of  ZrB2 
increases  with  temperature  (Fig.  11 -(b))  indicating  a  fast  redistribution  and  disposal  of  energy  in 
phonon  modes  in  ZrB2.  The  electronic  thermal  conductance  is  also  expected  to  rise  (shown  later). 
In  the  case  of  SiC,  the  non-availability  of  free  electrons  and  strong  covalent  bonds  associated 
with  possible  increase  in  vibration  amplitude  as  well  as  vibration  frequencies  due  to  increased 
temperatures  (shown  later)  should  lead  to  increase  in  phonon  thermal  conductance.  However,  as 
shown  the  phonon  thermal  conductance  reduces  with  increase  in  temperature.  Due  to  rigid  sp 3 
hybridized  bonds,  the  increased  thermal  energy  at  higher  temperatures  leads  to  increase  in 
phonon  frequencies.  Much  of  this  increase  goes  to  higher  frequency  optical  phonon  modes  that 
do  not  participate  in  thermal  conduction.  The  thermal  energy  in  acoustic  phonon  modes  that 
govern  thermal  conduction  reduces  leading  to  reduction  in  phonon  thermal  conductance  of  SiC. 
Since  SiC  does  not  have  free  electrons  its  electron  thermal  conductance  is  also  expected  to 
reduce  as  a  function  of  temperature  (shown  later).  Straining  in  both  cases  of  SiC  and  ZrB2  leads 
to  reduction  in  phonon  frequencies  relative  to  0%  strain  resulting  in  reduction  of  phonon  thermal 
conductance  with  increase  in  straining.  At  the  ZrB2-SiC  interface  the  thermal  conduction  is 
dominated  by  bonds  that  are  predominately  sp 3  hybridized  (only  Zr  related  bonds  have  metallic 
character)  leading  to  reduction  in  phonon  thermal  conductance  with  increase  in  temperature.  The 
value  of  the  interface  phononic  contribution  to  thermal  conductance  lies  close  to  the  ZrB2  value 
indicating  that  most  of  the  thermal  energy  at  the  interface  gets  redistributed  to  free  electrons  and 
ionic  cores  participating  in  thermal  conduction.  In  essence,  the  interface  has  attribute  of  SiC  with 
phononic  thermal  conductance  reducing  with  increase  in  temperature  while  the  thermal 
conductance  mechanism  is  that  of  ZrB2  with  thermal  energy  getting  redistributed  to  ionic  cores 
and  free  electrons. 

An  important  indicator  of  the  change  in  lattice  vibration  amplitude  and  frequencies  of  a 
material  is  its  thermal  diffusivity.  The  phononic  contribution  to  thermal  diffusivity  is  plotted  in 
Fig.  11 -(b)  at  0%  strain  value  as  a  function  of  temperature  for  all  three  crystals  examined.  As 
shown,  with  initial  increase  in  temperature  thermal  diffusivity  of  all  samples  follows  the  similar 
trend  as  that  of  the  phononic  thermal  conductance.  Later  in  the  case  of  ZrB2,  thermal  diffusivity 
starts  to  reduce  with  increase  in  temperature.  Thermal  diffusivity  is  an  indicator  of  how  fast  a 
material  can  achieve  thermal  equilibrium.  It  incorporates  thermal  conductivity  as  well  as  specific 
heat  capacity  of  a  material.  Faster  a  material  achieves  thermal  equilibrium  faster  it  is  able  to 
dispose  of  any  additional  heat  impingent  it.  Overall  trend  in  the  figure  is  of  reducing  thermal 
diffusivity  with  increase  in  temperature  as  well  as  with  increase  in  strain. 

The  observations  in  Fig.  11  and  related  arguments  are  intricately  linked  with  phononic 
vibration  frequencies  and  average  phonon  amplitude.  To  examine  the  issue  further,  average 
phonon  amplitude  of  all  systems  s  a  function  of  temperature  increase  and  strain  change  is 
examined  in  Fig.  12.  With  increase  in  temperature,  phonon  vibration  amplitude  increases  in  all 
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systems.  While  for  SiC  and  ZrB2,  the  increase  is  substantial  both  as  a  function  of  temperature 
and  strain,  the  interface  has  a  modest  increase.  Greater  thermal  energy  with  increase  in 
temperature  and  additional  strain  energy  due  to  increased  strain  energy  either  can  get  absorbed  in 
phonon  frequency  increase  or  in  phonon  vibration  amplitude  increase.  In  the  case  of  SiC  and 
ZrB2  the  additional  energy  changes  average  phonon  amplitude  significantly.  However,  due  to 
structural  mismatch  at  the  interface,  the  imposition  of  thermal  energy  or  strain  energy  does  not 
change  phonon  amplitude  at  the  interface  appreciably. 


Temperature,  [K]  Temperature,  [K] 

(a)  (b) 


Fig.  12  (a)  Ratio  of  average  phonon  amplitude  at  500  K  with  that  at  higher  temperatures  in  the 
case  of  ZrB2,  SiC,  and  ZrB2-SiC  interface  and  (b)  percentage  change  in  average  phonon 
amplitude  with  straining  from  0%  to  10%  as  a  function  of  temperature  in  the  case  of  ZrB2,  SiC, 
and  ZrB2-SiC  interface 

A  combination  of  Figs.  10  to  12  indicates  that  temperature  related  softening  of  examined 
materials  could  be  associated  with  increase  in  average  interatomic  distance  and  overall  reduction 
in  thermal  diffusivity.  With  material  softening  at  higher  temperatures  and  the  corresponding 
increase  in  interatomic  distance  the  average  phonon  amplitude  increases  that  leads  to  higher 
thermal  energy  absorption  capability.  Straining,  invariably,  led  to  change  in  phonon  amplitude  as 
well  as  phonon  density  of  states.  Both  factors  combined  together  to  further  reduce  the  thermal 
diffusivity  in  materials  as  a  function  of  temperature.  These  observations  also  indicate  that  by 
phase  mixing  it  is  possible  to  tune  thermal  conductance  as  well  as  rate  of  conductive  heat 
transport  in  a  material. 

Besides  phonon,  electrons  have  been  considered  to  be  of  importance  in  thermal  conduction 
especially  at  high  temperatures,  [59].  The  electron  contribution  to  thermal  conduction  in  the  form 
of  electron  thermal  conductance  was  calculated  using  the  NEGF  formulation  described  before, 
Fig.  13.  The  trend  of  the  electron  thermal  conductance  variation  as  a  function  of  temperature  and 
strain  remains  the  same  as  observed  in  the  case  of  the  phonon  thermal  conductance  in  the  case  of 
SiC  and  ZrB2.  However,  interface  now  shows  increase  in  electron  thermal  conductance  with 
increase  in  temperature.  The  variation  with  temperature  increase  shown  by  the  interface  electron 
thermal  conductance  is  also  significantly  higher  in  comparison  to  the  one  shown  by  interface 
phonon  thermal  conductance.  Electron  thermal  diffusivity  also  shows  a  similar  trend  as  phonon 
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thermal  diffusivity  (not  shown).  With  softening,  the  electron  thermal  conductance  shows 
appreciable  increase  in  the  case  of  interface  and  ZrB2.  A  comparison  of  thermal  diffusivity  of 
electrons  and  phonons  separately  at  high  temperatures  indicates  that  thermal  diffusivity 
contribution  of  electrons  increases  at  high  temperatures  to  as  much  as  20%  in  comparison  of 
about  10%  at  room  temperature.  These  observations  are  in  line  with  earlier  discussion  regarding 
availability  of  free  electrons  in  the  case  of  ZrB2  and  interface  and  non-availability  of  such 
electrons  in  the  case  of  SiC.  With  increase  in  temperature  electrons  get  increasingly  mobile  and 
with  reduced  bad  gap  start  to  participate  significantly  in  the  thermal  conduction.  Interface 
formation  between  SiC  and  ZrB2  leads  to  stronger  material  with  reduced  thermal  diffusivity  but 
higher  electron  thermal  conductance. 
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Figure  13  (a)  Electron  thermal  conductance  of  single  crystalline  SiC,  single  crystalline  ZrB2, 
and  SiC-ZrB2  interface  as  a  function  of  strain  and  temperature  and  (b)  ratio  of  average  electron 
energy  at  500  K  with  that  at  higher  temperatures  in  the  case  of  ZrB2,  SiC,  and  ZrB2-SiC  interface 

Contribution  of  ZrB2  to  the  interface  electron  density  of  states  is  examined  in  Fig.  14  using 
Mulliken  analysis.  The  figure  shows  change  in  J-orbital  electron  density  of  states  as  a  function  of 
temperature  at  0%  strain.  Appreciable  change  in  electron  density  is  found  to  occur  as  a  function 
of  temperature  in  the  case  of  interface  in  comparison  to  that  in  the  case  of  ZrB2.  With  ZrB2 
structural  setup  any  influence  of  atomic  restructuring  occurs  by  shift  in  the  lighter  B  atoms.  In 
the  case  of  interface,  besides  Zr-B  bonds  (combined  sp2 -sp 3  hybridized  with  spd  contributions), 
Si-B  (sp3  hybridized)  bonds  also  participate  in  shift  of  B  atoms.  This  behavior  of  lighter  atoms 
taking  more  shifts  under  the  influence  of  external  parameters  to  adjust  for  the  changes  in  the 
surrounding  may  be  one  factor  that  can  contribute  to  increasing  average  electron  energy  trend 
with  respect  to  temperature  in  the  case  of  interface.  Another  aspect  is  that  interface  d-orbital 
electron  density  of  states  is  distributed  over  negative  as  well  as  positive  energy  levels.  The  ZrB2 
electron  density  of  states  is  distributed  over  positive  energy  levels  predominantly.  Clearly,  the 
delocalized  electrons  observed  in  electron  density  contour  at  the  SiC  side  of  interface  in  Fig.  15, 
have  significant  influence  of  valence  J-orbital  in  the  ZrB2  phase  of  the  interface.  Such 
delocalization  imparts  metastability  to  interface  and  ultimately  leads  to  increase  in  average 
electron  energy  as  a  function  of  increase  in  temperature  as  well  as  increase  in  electron  thermal 
conductance  of  interface. 
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(a)  (b) 

Fig.  14  Mulliken  analysis  based  d-orbital  electron  density  of  states  comparison  at  500  K  and 
2500  K  and  at  0%  strain  in  the  case  of  (a)  ZrB2  single  crystal  and  (b)  ZrB2-SiC  interface 


Fig.  15  Electron  density  contour  for  ZrB2-SiC  interface  at  (a)  300  K  and  0%  strain,  (b)  300  K 
and  10%  tensile  strain  (c)  at  2500  K  and  0%  strain  and  (d)  at  2500  K  and  10%  strain. 
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§2.1.1  SUMMARY 

Fundamentally,  tunable  thermal  shock  performance  of  a  high  temperature  material  demands 
high  thermal  diffusivity  combined  with  high  thermal  conductivity.  A  first  step  in  this  direction  is 
to  be  able  to  examine  an  interface  strength  at  high  temperatures  and  its  correlation  with  thermal 
conductivity  and  thermal  diffusivity.  Since,  high  temperature  operation  involves  possible 
complex  nanostructural  phase  changes  in  an  interface,  a  quantum  mechanical  approach  to 
analyze  such  correlation  is  necessary.  The  present  work  reported  such  analyses.  The  most  stable 
<0001>-<111>  ZrB2-SiC  interface  is  chosen  for  analyses  as  a  function  of  strain  and  temperature 
for  correlation  between  the  thermal  and  mechanical  properties. 

It  is  widely  known  that  with  temperature  increase,  ratio  of  individual  contribution  of  phonons 
and  electrons  to  overall  thermal  conductance  changes,  but  very  little  is  known  to  what  extent. 
The  present  work  sheds  some  light  in  this  regard  for  SiC  and  ZrB2  single  crystals  as  well  as  for 
the  SiC-ZrB2  interface.  Earlier,  Marshall  group  [59]  have  analyzed  thermal  conduction  at  high 
temperatures  in  a  variety  of  ZrB2-SiC  composites.  However,  their  findings  results  are  at  the 
composite  macro  level.  Findings  in  this  work  regarding  interface  and  single  crystal  thermal- 
mechanical  correlations  can  significantly  help  in  interpreting  their  results,  especially  with  regard 
to  the  role  played  by  electron  thermal  conductance.  The  findings  in  this  work  for  the  first  time 
analyze  thermal-mechanical  correlations  in  single  crystalline  phases  and  compare  those  with  the 
behavior  of  an  interface. 

Analyses  indicate  that  there  is  a  direct  correlation  between  strength  reduction  with  increase  in 
temperature  and  the  corresponding  proportional  increase  in  average  interatomic  distance  in  the 
case  of  all  three  samples.  The  phonon  thermal  diffusivity  shows  a  direct  and  proportionate 
correlation  with  temperature  increase  in  the  case  of  ZrB2  and  an  inverse  correlation  in  the  case  of 
SiC  and  interface,  indicating  a  possibility  that  in  these  materials  temperature  increase  may  be 
associated  with  the  changes  in  the  ability  to  diffuse  heat.  Thermal  conductance  is  found  to  have  a 
strong  correlation  to  phonon  density  of  states.  The  thermal  diffusivity,  however,  is  found  to  be 
correlated  with  both  phonon  amplitude  and  phonon  frequencies.  Fundamentally,  increased 
phonon  amplitude  correlates  to  increase  in  specific  heat  and  therefore  reduction  in  thermal 
diffusivity.  Effect  of  strain  is  always  to  reduce  phonon  density  of  states  and  increase  the  phonon 
vibration  amplitude.  Reduction  in  the  density  of  states  leads  to  reduced  thermal  conductance  and 
increase  in  phonon  vibration  amplitude  leads  to  increased  specific  heat  and  reduced  thermal 
diffusivity.  Analyses  of  electron  thermal  conductance  and  electron  energy  change  with  respect  to 
temperature  indicate  that  electron  thermal  contributions  increase  appreciably  with  temperature. 
While  electrons  get  delocalized  with  increase  in  temperature  and  strain,  their  average  energy 
occupation  levels  increase  leading  to  reduction  in  their  average  energy  in  the  case  of  single 
crystals.  However,  the  delocalization  does  not  lead  to  reduction  in  electron  energy  levels  at  the 
examined  interface  indicating  a  metastable  state.  Examination  of  change  in  thermal  properties  at 
different  mechanical  strain  levels  reveals  that  the  mechanisms  of  strength  and  thermal  property 
change  with  increase  in  temperature  may  be  similar  to  the  mechanisms  responsible  for  property 
change  with  change  in  applied  strain. 

§2.2  FINDINGS  FROM  CLASSICAL  MOLECULAR  SIMULATIONS 

.  Since  ab-initio  approach  is  limited  to  small  system  sizes  because  of  significant 
computational  demand,  such  simulations  are  being  coupled  to  approximate  classical  molecular 
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simulations  to  scale  up  the  simulation  predictions  to  realistic  microstructure  size.  Accordingly, 
the  development  of  ab-initio  approach  was  performed  together  with  analyses  of  interfacial  heat 
transfer  as  a  function  of  molecular  scale  morphology.  A  widely  studied  material  system  (Si-Ge) 
was  analyzed  in  order  to  standardize  the  approach  in  light  of  available  data.  In  the  following 
salient  details  of  the  simulations  are  presented. 

§2.2.1  THERMAL  CONDUCTION  ANALYSES  IN  SUPERLATTICES  AS  A  FUNCTION 
OF  STRAIN 

Superlattices  are  nanoscale  engineered  material  system  in  which  the  thermal  conduction 
properties  could  be  tailored  for  applications  such  as  high  figure  of  merit  (ZT)  thermoelectric, 
microelectronics,  and  optoelectronics  devices  etc.  [60]  Widely  researched  superlattices  such  as 
Si/Sii_xGex,  [61-63]  GaAs/AlAs,  [64]  and  Bi2Te3/Sb2Te3,  [65]  have  been  observed  to  have 
thermal  conductivity  values  much  lower  than  the  bulk  systems  of  similar  composition.  Factors 
that  could  be  adjusted  for  tailoring  the  thermal  conductivity  of  superlattices  include  the 
monolayer  film  thickness,  periodicity,  heat  flow  direction,  straining,  and  temperature  of 
operation.  Different  groups,  [61,  62],  have  worked  on  experimentally  analyzing  thermal  transport 
in  Si-Ge  superlattice  thin  film  structures  as  a  function  of  the  number  of  superlattice  period  and 
thickness.  Chen  and  coworkers  analyzed  thermal  conductivity  in  superlattice  systems  using 
numerical  solutions  of  Boltzmann  transport  equations  (BTE)  and  molecular  dynamics  (MD). 
[66]~[67]  Chen  and  co-workers  [68]  performed  MD  simulations  studying  the  dependence  of 
thermal  conductivity  on  the  period  length.  NEMD  analyses  of  Si-Ge  superlattices  with 
unspecified  number  of  periods  using  Stillinger-Weber  potential  have  been  performed  to  study  the 
dependence  of  thermal  conductivity  on  the  monolayer  film  thickness.  [69,  70] 

Thus,  analyses  so  far  lack  a  combined  account  of  factors  such  as  tensile  vs.  compressive 
straining,  change  in  number  of  periods,  change  in  temperature,  and  change  in  the  monolayer  film 
thickness  for  a  single  superlattice  material  system.  Recently,  Zhou  and  coworkers  [71]  did  an 
analysis  of  the  effect  of  heat  flow  direction  on  thermal  conductivity  of  composite  metals.  They 
found  that  the  thermal  properties  change  with  the  change  in  heat  flow  direction.  The  present 
work  attempts  to  analyze  the  effect  of  tensile  vs.  compressive  straining,  change  in  number  of 
periods,  change  in  temperature,  and  change  in  the  monolayer  film  thickness  on  the  thermal 
conductivity  of  Si-Ge  superlattices. 


(5Six5Ge)3:  Three  Periods  of  5  nm  Si  &  5  nm  Ge 
Si  Ge 


Fig.  16:  Terminology  for  layered  structures 
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NON-EQUILIBRIUM  MOLECULAR  DYNAMICS 

Figure  16  shows  the  simulation  setup  and  the  terminology  for  the  superlattices  analyzed  in 
the  presented  research  using  NEMD.  The  simulation  supercell  is  bounded  by  fixed  atomic  layers 

O 

of  length  equivalent  to  25  A,  i.e.  up  to  9  atomic  layers  on  either  side,  Fig.  16.  In  order  to  simulate 
a  thin  film  structure,  PBCs  are  imposed  in  the  directions  transverse  to  the  length  of  supercells. 
All  layers  have  <100>  orientation  along  the  length  of  supercells.  We  chose  4x4  unit  cells  in 
transverse  direction  cross-sectional  area,  [72],  As  shown,  a  (5six5ce)3  structure  represents  an 
interfacial  supercell  of  5  nm  thick  Si  thin  film  with  5  nm  thick  Ge  thin  film.  Subscript  3  denotes 
that  the  supercell  consists  of  three  periods.  In  the  presented  research,  supercells  with  four 
different  thicknesses  (2.5nm,  5nm,  7.5nm,  lOnm)  and  three  different  periods  (1,  2,  3)  are 
analyzed  at  three  different  temperatures  (400  K,  600  K,  and  800  K).  We  can  calculate  thermal 
conductivity,  k,  as 


k  = 


Ag/i<»  + 

tsA  |v  T 


(8) 


Here, vr  =  (a til)  where  at  is  the  temperature  difference  between  the  hot  and  cold  reservoirs, 

L  is  the  length  of  simulation  cell,  A  is  the  cross-sectional  area  of  the  simulation  supercell  and  ts  is 
the  simulation  time.  To  impose  a  constant  heat  flux,  the  hot  and  cold  reservoirs  need  to  be  kept  at 
a  constant  temperature.  The  energy  supplied  to  the  hot  and  cold  reservoirs  to  maintain  a  constant 
temperature  is  based  on  momentum  conservation  scheme  [73]  and  can  be  calculated  from 


A  E 


2 


)  • 


(9) 


Here,  v"’“  is  the  old  atomic  velocity,  v‘-new  is  the  new  atomic  velocity  after  scaling.  In  addition 
to  computing  k,  thermal  boundary  resistance  (TBR)  across  an  interface  i,  r‘bd  can  be  calculated 
as 


rIbd  = 

iE>9> 

A  ts 

Here,  A  rint  represents  the  temperature  drop  measured  at  each  interface  and  J  gives  the  heat  flux. 

To  study  the  effect  of  strain  on  the  thermal  conductivity,  we  varied  the  strain  by  stretching  the 
simulation  cell  from  10%  to  -10%  with  step  size  of  2%  and  ran  MD  simulation  at  each  strain 
level  to  calculate  the  thermal  conductivity. 

Simulation  Setup 

The  inter- atomic  interactions  for  Si-Ge  systems  are  described  by  the  Tersoff  bond-order 
potential,  [74],  During  NEMD  simulations  to  compute  k  at  a  temperature  T,  the  superlattice 
system  is  equilibrated  for  200  ps  with  a  time  step  of  1  femtosecond  (fs)  in  microcanonical  (NVE) 


(10) 

(11) 
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ensemble  at  temperature  T.  After  equilibration,  a  temperature  gradient  is  established  by  imposing 
Thot=T+  30  K  and  TC(M=T-?>()  K  in  the  hot  and  the  cold  reservoir  respectively  using  momentum 
conservation  scheme  [73],  followed  by  further  equilibration  of  the  computational  supercell  for 
500  ps.  During  this  equilibration  procedure,  at  each  time  step  the  values  of  k  are  calculated  using 
Eq.  (8).  Calculations  showed  that  the  heat  flux  imposed  on  the  superlattices  by  fixing  up  the  hot 
and  cold  reservoir  temperatures  took  approximately  from  200  to  300  ps  to  get  stabilized.  For 
calculating  temperature  profile  along  the  supercell  length,  each  supercell  was  divided  into  thin 
slabs  of  length  a  little  larger  than  both  Si  and  Ge  lattice  constant.  Once  the  values  of  k  converge, 
the  temperature  profile  along  the  length  of  the  supercells  is  obtained  by  calculating  average 
temperature  of  each  slab  based  on  the  total  kinetic  energy  of  all  atoms  in  the  slab  averaged  over 
100  ps. 


Fig.  17:  Temperature  profile  for  (a)  (10SiXl0ge)i  and  (10gexl0Si)i,  and  (b)  (10SiXl0ge)2  and 
(10gexl0si)2  simulation  systems 

RESULTS 

Figure  17  displays  the  temperature  profile  obtained  along  the  length  of  the  (10siXl0Ge)x=i,  2 
and  (lOceXlOsiW,  2  simulation  systems  plotted  for  increasing  period  values.  Fig.  2(a)  shows 
temperature  profile  for  (10siXl0Ge)i  and  ( 1  OceX  1  Osi )  1  superlattices  after  convergence  in  k  values 
corresponding  to  F=400  K.  As  shown  in  the  curves  in  Fig.  17(a),  a  steep  drop  in  temperature  is 
observed  at  the  interface  of  Si-Ge  at  the  position  corresponding  to  10  nm.  Such  drop  in  the 
temperature  is  attributed  to  the  thermal  boundary  resistance  offered  by  the  interface  and  is 
observed  to  be  different  for  the  superlattice  with  different  directions  of  heat  current.  This 
directional  dependence  is  discussed  later  in  section.  Fig.  17  (b)  shows  the  temperature  drop 
across  interfaces  of  (10siXl0Ge)2and  (lOccXlOsiE  which  also  clearly  shows  the  effect  of  interface 
boundary  resistance.  Another  important  aspect  of  interfacial  conduction  observed  in  Fig.  17  is 
the  non-linear  behavior  of  interfaces  in  offering  resistance  to  heat  flow  as  the  number  of 
interfaces  increase.  Calculations  done  on  the  basis  of  temperature  drop  recorded  at  each  interface 
leads  to  the  observation  that  the  total  TBR  does  not  increase  linearly  with  the  increase  in  the 
number  of  interfaces,  as  suggested  by  Ref.  [75], 
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Thermal  Conductivity  as  a  Function  of  Superlattice  Period,  Film  Thickness,  Temperature 
and  Strain 

We  examined  the  variation  of  thermal  conductivity  of  (5six5Ge)i,2,  and  3  and  (10sixl0Ge)i,2,  and  3 
superlattices  with  increase  in  number  of  periods  at  temperatures  400K,  600K  and  800K.  Thermal 
conductivity  increases  with  increase  in  monolayer  thickness  and  with  increase  in  number  of 
periods.  The  rate  of  change  of  thermal  conductivity  with  increase  in  the  number  of  periods  is 
found  to  be  higher  for  thicker  monolayer  films.  This  is  because,  there  are  two  factors  competing 
with  each  other:  (1)  number  of  interfaces  increase  leading  to  higher  cumulative  rbd  ,  and  (2) 

length  of  the  superlattice  system  increase  leading  to  a  drop  in  the  overall  temperature  gradient. 
For  higher  superlattice  period  thickness,  i.e.  with  thicker  monolayer  films,  the  latter  dominates 
over  the  former  by  a  significant  proportion.  Accordingly,  there  is  a  steeper  increase  in  the 
thermal  conductivity  value  for  higher  period  thickness. 


Fig.  18:  Thermal  conductivity  as  a  function  of  strain  at  different  temperatures  for  (5six5ge)  3, 
(10Sixl0ge)  3  structures 

Effect  of  Straining  on  Thermal  Conductivity  as  a  Function  of  Temperature 

Figure  18  displays  the  thermal  conductivity  of  (5six5ce)3  and  (lOsixlOceri  systems  as  a 
function  of  temperature  at  different  strain  levels.  Both  compressive  and  tensile  strain  causes  the 
thermal  conductivity  of  superlattices  to  decrease.  This  trend  is  different  from  the  earlier  reported 
analyses  in  literature  [75-77]  regarding  the  effect  of  strain  on  the  thermal  conductivity.  We 
conjecture  that  this  difference  is  attributable  to  simulations  being  performed  near  or  above  Debye 
temperature  values  in  our  case.  With  the  increase  in  the  number  of  periods,  the  decrease  in 
thermal  conductivity  is  higher  for  higher  strain  and  this  trend  is  observed  at  all  three 
temperatures:  400K,  600K  and  800K.  It  is  observed  that  the  straining  has  insignificant  effect  on 
the  thermal  conductivity  of  superlattices  with  5  nm  monolayer  thickness,  as  we  approach  thin 
film  limit.  Similarly,  when  tensile  strain  is  applied,  atomic  distance  increases,  phonon  relaxation 
time  increases  and  structure  stiffness  decreases.  But,  an  increase  in  the  period  thickness  and 
periodicity  causes  the  nucleation  of  structural  defects,  [62]  (identified  as  deviation  from  perfect 
single  crystal  structure).  The  formation  of  structural  defects  is  more  pronounced  at  high 
temperatures,  which  is  experienced  in  this  study.  The  scattering  at  the  defects  and  dislocations 
provide  additional  resistance  to  the  flow  of  heat  current  through  the  material  system,  which  leads 
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to  a  reduction  in  the  thermal  conductivity  with  an  increase  in  tensile  strain,  [75],  It  is  also  clear 
from  the  figures  that  the  effect  of  straining  is  more  pronounced  at  larger  period  thickness,  which 
is  again  due  to  chances  of  dislocation  nucleation  being  higher  for  higher  period  thickness.  This 
decrease  of  thermal  conductivity  due  to  strain  is  seen  to  increase  with  an  increase  in  the  number 
of  periods,  at  all  period  thicknesses. 

Effect  of  Heat  Flow  Direction  Reversal  on  Thermal  Conductivity  Values 

Figure  19  shows  a  comparison  of  thermal  conductivity  as  a  function  of  temperature  and  heat 
flow  direction  for  (7.5siX7.5Ge)3  and  (7.5Gex7.5si)3  superlattices,  and  for  (10siXl0Ge)3  and 
(10Gexl0si)3  superlattices.  As  shown,  with  an  increase  in  the  period  thickness,  the  effect  of  the 
reversal  in  the  heat  flow  direction  becomes  significant.  This  behavior  can  be  attributed  to  the 
change  in  the  frequency  of  the  heat  carrying  phonons  with  the  change  in  the  heat  flow  direction. 
For  Si— >Ge  system,  heat  transfer  characteristics  are  dominated  by  the  phonons  in  Si  monolayer 
whereas  in  Ge— >Si,  the  heat  transfer  is  determined  by  the  phonons  in  Ge  monolayer.  Owing  to  a 
large  atomic  mass  difference  between  Si  and  Ge,  which  reflects  both  in  terms  of  large  acoustic 
mismatch  across  Si-Ge  layer  and  also  in  the  difference  in  the  frequency  of  Si  and  Ge  phonons, 
the  thermal  resistance  offered  by  an  interface  to  a  particular  type  of  phonon  varies.  This 
difference  grows  with  an  increase  in  the  number  of  periods  and  results  in  a  larger  drop  in  thermal 
conductivity  value  for  Ge— >Si  system,  when  compared  with  Si— >Ge  system. 


Fig.  19:  Dependence  of  thermal  conductivity  on  heat  flow  direction  as  a  function  of  temperature 
for  (a)  (7.5six7.5ge)4  and  (7.5gex7.5si)4  superlattices,  and  (b)  (10siXl0ge)4  and  (10gexl0si)4 
superlattices 

§2.2.2  THERMAL  CONDUCTION  IN  BIOMIMETIC  COMPOSITES  AS  A  FUNCTION 
OF  STRAIN 

The  MD  based  thermal  conductivity  calculations  in  the  case  of  biomimetic  nanocomposites  are 
based  on  the  same  framework  used  in  the  case  of  superlattices.  Biomimetic  nanocomposites  with 
nanostructures  in  the  form  of  nanoparticles  and  nanowires  embedded  in  a  host  matrix  material 
differ  from  conventional  composite  materials  in  their  thermal  and  mechanical  properties  due  to 
exceptionally  high  surface  to  volume  ratio  of  the  nanostructures  acting  as  reinforcing  phase. 
Atomistic  analyses  of  the  biomimetic  nanocomposite  thermal  behavior  are  limited.  Jeng  and 
coworkers  [78]  have  used  Monte  Carlo  simulations  to  study  the  phonon  transport  and  thermal 
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conductivity  reduction  in  the  biomimetic  nanocomposites  with  Si  nanoparticle  embedded  in  Ge 
host.  Similar  Monte  Carlo  analyses  for  the  nanocomposites  with  Si  tubular  nanowires  in  Ge  host 
have  been  performed  by  Yang  and  group  [79].  In  a  recent  work,  Huang  and  group  [80]  have 
used  molecular  dynamics  (MD)  simulations  to  analyze  similar  systems  as  analyzed  by  Yang  and 
group  [79].  Yang  and  group  [81]  have  theoretically  studied  phonon  thermal  conductivity  of 
periodic  two-dimensional  nanocomposites  with  nanowires  embedded  in  a  host  semiconductor 
material  using  phonon  Boltzmann  equation.  They  have  concluded  that  the  thermal  conductivity 
of  nanocomposites  is  always  higher  than  that  of  a  superlattice  with  the  same  characteristic 
thickness.  The  difference  reduces  as  the  dimensions  of  the  two  types  of  materials  are  reduced.  A 
single  conclusion  common  to  all  analyses  has  been  that  the  reason  behind  biomimetic 
nanocomposites  displaying  low  thermal  conductivities  is  high  percentage  of  atoms  in  the 
interfaces  leading  to  significant  phonon  scattering.  Interface  atom  fraction  in  the  superlattices  is 
usually  lower  than  that  in  the  biomimetic  nanocomposites.  Since  structurally  the  superlattices 
and  the  biomimetic  nanocomposites  are  significantly  different,  their  response  to  externally 
imposed  mechanical  straining  should  also  be  different.  By  a  judicious  combination  of  applied 
straining  and  morphology  it  should  be  possible  to  tailor  the  thermal  conductivity  of  both  types  of 
materials.  With  this  view,  the  present  investigation  focuses  on  understanding  thermal  behavior 
of  a  set  of  Si-Ge  biomimetic  nanocomposites  using  non-equilibrium  molecular  dynamics 
(NEMD)  simulations  at  three  different  temperatures  (400  K,  600  K,  and  800  K)  and  at  strain 
levels  varying  between  -10%  and  10% 


Fig.  20  Simulation  setup  and  terminology  for  the  nanocomposite  structures  analyzed  for  thermal 
conduction.  The  staggered  arrangement  of  Ge  blocks  inside  Si  host  matrix  is  displayed  along 
with  various  parameters  considered  in  this  study 

Figure  20  displays  the  staggered  arrangement  of  Ge  blocks  inside  Si  host  matrix,  highlighting 
various  parameters  that  were  varied  to  study  their  effect  on  the  overall  thermal  conductivity  of 
the  biomimetic  nanocomposites.  The  nomenclature  (SiiooGe°6)4oo  shown  in  the  Fig.  20  helps  in 
identifying  each  biomimetic  nanocomposite  analyzed.  The  subscript  100  on  Si  represents  the 
thickness  of  the  nanocomposite  with  Ge  cubic  blocks  of  size  6  A  as  indicated  by  the  Ge  subscript 
of  6.  The  superscript  0  to  Ge  represents  the  Y-shift  (extent  to  which  a  Ge  particle  is  shifted  in  the 
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y-direction  with  respect  to  an  adjacent  Ge  particle).  In  the  present  work,  nanocomposites  with 
three  different  thicknesses:  lOnm,  20nm  and  30nm  are  analyzed.  For  each  thickness, 
nanocomposites  with  three  different  Ge  block  sizes:  6  A,  9  A,  and  12  A  are  generated.  Four 

o  o  o  o 

values  of  Y-shift:  0  A,  2  A,  4  A,  and  6  A  are  used  for  each  nanocomposite  thickness  and  for  each 
Ge  block  size.  It  is  made  sure  that  all  nanocomposites  have  equal  number  of  Ge  blocks  at  a 
cross-sectional  view  in  the  x-y  plane.  As  shown  in  Fig.  20,  the  simulation  supercell  is  bounded 
by  fixed  atomic  layers  of  thickness  equivalent  to  25A,  i.e.  up  to  9  atomic  layers  on  either  side. 
Further  increase  in  the  thickness  of  fixed  atomic  layers  did  not  change  the  presented  results. 
Periodic  boundary  conditions  (PBCs)  are  imposed  in  the  directions  transverse  to  the  thickness  of 
composite  supercells  to  be  able  to  simulate  a  thin  film  structure.  Both,  Si  and  Ge  have  <100> 
orientation  along  the  thickness.  A  nanocomposite  configuration  may  affect  the  lattice  constant 
values  of  both  Si  and  Ge.  Based  on  the  analyses  by  Volz  and  coworkers  [82],  lattice  parameters 
for  Si,  asi=5.43  A  and  for  Ge,  ace=5.657  A  are  used. 

Effect  of  Straining  as  a  Function  of  Temperature 


(a) 


(b) 


(c) 


(d) 


Fig.  21  Comparisons  of  the  thermal  conductivity  values  of  nanocomposites  and  superlattices  of 
equal  thickness  as  a  function  of  strain.  Plots  are  (a)  for  Si|ooGe4y  and  (5six5Ge)i  under  tension,  (b) 
for  SiiooGe49and  (5six5ce)i  under  compression,  (c)  for  Si2ooGe49  and  (10sixl0ce)i  under  tension, 
and  (d)  for  Si2ooGe4g  and  (lOsixlOGe)i  under  compression  at  400  K,  600  K  and  800  K. 

Figure  21  displays  the  thermal  conductivity  variation  of  SiiooGe49  and  Si2ooGe49 
nanocomposites  as  a  function  of  tensile  and  compressive  straining  at  three  different 
temperatures:  400K,  600K  and  800K.  Variation  of  the  thermal  conductivity  of  the  superlattices: 
(5six5Ge)i  and  (10sixl0Ge)i,  [83],  with  comparable  thickness  has  been  plotted  (dashed  lines)  along 


35 


the  nanocomposite  thermal  conductivity  values  in  order  to  offer  a  comparison.  As  shown,  a 
linear  fit  to  the  data  on  thermal  conductivity  for  the  nanocomposites  as  well  as  the  superlattices 
can  be  obtained.  It  is  observed  that  the  straining  has  insignificant  effect  on  the  thermal 
conductivity  of  the  nanocomposites  with  10  nm  thickness  and  slightly  affects  the  thermal 
conductivity  of  the  nanocomposites  with  20  nm  thickness.  Overall,  however,  the  thermal 
conductivity  shows  a  stronger  dependence  on  strain  in  the  case  of  superlattices  when  compared 
to  the  nanocomposites,  [83],  The  superlattices  with  comparable  thickness  show  an  increase  in 
thermal  conductivity  values  as  a  function  of  tensile  strain  and  decrease  in  the  values  as  a  function 
of  compressive  strain.  This  difference  of  the  effect  of  straining  on  the  thermal  conductivity  of 
nanocomposites  when  compared  with  superlattices  can  be  explained  from  the  phonon  spectral 
density  plots.  In  superlattice  phonon  spectrum,  (not  shown  here,  can  be  observed  in  our  Journal 
of  Physics  Publication  pointed  out  at  first  page),  we  have  sharply  defined  peaks  corresponding  to 
certain  frequencies  which  dominate  the  heat  transfer.  Correspondingly,  only  a  limited 
frequencies  are  the  most  significant  contributors  to  the  heat  transfer.  On  the  contrary,  in  the 
nanocomposites  there  are  no  sharply  defined  peaks  in  the  phonon  spectrum.  Accordingly,  there  is 
a  wide  spectrum  of  phonon  frequencies  with  equivalent  spectral  densities,  which  carry  the  heat 
current  across.  Therefore,  the  reduction  or  increase  in  the  interatomic  distance  caused  by 
straining  and  the  ensuing  gain  or  loss  in  the  stiffness,  respectively,  can  affect  superlattice  phonon 
spectrum  more  owing  to  the  dominance  of  a  limited  frequencies.  In  the  nanocomposites,  the 
effect  is  subdued,  because  the  heat  transfer  is  distributed  across  a  much  larger  spectrum  of 
wavelengths,  Fig.  22.  This  difference  of  heat  carrier  phonons  in  the  nanocomposites,  when 
compared  to  those  in  superlattices  leads  to  the  observed  difference  in  the  variations  of  thermal 
conductivity  shown  in  Fig.  6.  To  further  explain  the  variation  of  thermal  conductivity  with 
straining,  we  can  write  k  as,  [75,  77], 

k  =  —C  v2t  Or  k  =  —  C  r.  2  —  t  .  (12) 

0  Ph  0  1  J  Ph  V  7 

3  3  m 

Here  r  ph  gives  phonon-phonon  interaction  or  phonon  relaxation  time,  g  represents  the 
stiffness  constant  of  the  atomic  structure,  and  m  gives  the  mass  of  atoms.  In  the  above  equation, 
the  relations  i  =  vrph  and  v  =  r.Jg  Jm  are  used.  Here,  r  represents  the  average  interatomic 
distance.  It  can  be  shown  that  an  increase  in  the  compression  causes  a  decrease  in  r 
equivalently  for  both  superlattices  and  nanocomposites. 


Figure  22  A  comparison  of  the  phonon  spectral  density  at  different  strain  levels  for  SiiooGe49 
nanocomposite  at  400  K,  under  (a)  tension,  and  (b)  compression 
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As  shown  in  Fig.  23,  the  radial  distribution  functions  (RDFs)  for  the  nanocomposites  are  not 
significantly  affected  by  straining.  Therefore,  tension  or  compression  affects  superlattice  thermal 
conductivity  in  a  more  significant  manner.  It  is  also  clear  from  the  Fig.  21  that  the  effect  of 
straining  is  more  pronounced  at  larger  period  thickness  for  superlattices  and  higher  thickness 
nanocomposites,  which  can  be  attributed  to  chances  of  defect  nucleation  being  higher  for  higher 
period  thickness.  In  addition,  for  higher  thickness  nanocomposites,  the  interfacial  atom  fraction 
reduced  when  compared  to  that  in  smaller  thickness  nanocomposites  (e.g.  10  nm  here).  This 
directly  leads  to  higher  effect  of  straining  in  higher  thickness  nanocomposites. 


(a) 


(b) 


(c) 


Figure  23  A  comparison  of  the  (a)  Si-Si,  (b)  Ge-Ge,  and  (c)  Si-Ge  normalized  radial  distribution 
functions  (RDF)  at  different  strain  levels  for  SiiooGe4g  nanocomposite  at  400  K 

The  reduction  in  the  interatomic  distance  due  to  compression  or  increase  in  the  interatomic 
distance  due  to  tension  significantly  affects  phonon  relaxation  time,  t .  Such  changes  in  the 

interatomic  distances  also  lead  to  the  corresponding  changes  in  the  lattice  stiffness.  Since, 
nanocomposites  have  a  significantly  high  interface  atom  fraction,  the  extent  of  r..  reduction  or 

increase  with  tension  or  compression,  respectively,  is  limited  when  compared  to  that  in 
superlattices,  Fig.  23. 
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